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The competing orders in the particle-particle (P-P) channel and the particle-hole (P-H) chan- 
nel have been proposed separately to explain the pseudogap physics in cuprates. By solving the 
Bogoliubov-deGennes equation self-consistently, we show that there is a general complementary 
connection between the d-wave checkerboard order (DWCB) in the particle-hole (P-H) channel and 
the pair density wave order (PDW) in the particle-particle (P-P) channel. A small pair density 
localization generates DWCB and PDW orders simultaneously. The result suggests that suppress- 
ing superconductivity locally or globally through phase fluctuation should induce both orders in 
underdoped cuprates. The presence of both DWCB and PDW orders with 4a x 4a periodicity can 
explain the checkerboard modulation observed in FT-STS from STM and the puzzling dichotomy 
between the nodal and antinodal regions as well as the characteristic features such as non-dispersive 
Fermi arc in the pseudogap state. 



PACS numbers: 74.25.Jb, 74.25. Dw, 74.72.-h 



I. INTRODUCTION 

An important characteristic of strongly correlated elec- 
tron systems is the existence of different instabilities that 
lead to many competing orders, fn high temperature su- 
perconductors, besides the superconducting phase, many 
competing orders, such as, spin density wave (SDW)i£, 
d-density wave (DDW) 3,4 , pair density wave (PDW) 5,6 , 
stripe 7 -, and so on, have been proposed to explain vari- 
ous experimental observations. Those competing orders 
can be generally classified into two categories, the orders 
in particle-particle (P-P) channel and the orders in the 
particle-hole (P-H) channel. So far, most of theoretical 
works in cuprates have focused on the effect of individual 
competing orders. However, the orders in two channels 
are not completely independent of each other. In some 
cases, they must be correlated. In this paper, we detailed 
study one of these examples, the extended s-wave PDW 
order in the P-P channel and the <i-wave checkerboard 
density (DWCB) order in the P-H channel. 

The motivation of this study mainly comes from 
the recent experiments of Scanning Tunneling Mi- 
croscopy(STM) in cuprates. These experiments have 
revealed surprising yet important electronic structures 
in the high temperature superconductors. The Fourier 
transform scanning tunneling spectroscopies (FT-STS) 
from STM have captured two different general features in 
both momentum and energy space o 8 ' 9 ' 10 : 11 ] 12 ! 13 ] 14 : 15 : 16 . 
One feature is dispersive peaks in FT-STS 11 ' 12 , inter- 
preted as interference patterns caused by elastic scat- 
tering of quasiparticles from impurities 4 -. The other 
is non-dispersive peaks, a checkerboard modulation ob- 
served in various different materials and circumstances. 
The checkerboard structure was first discovered locally in 
BSCCO near a vortex core 8 -^. Then, it was found to be 
a characteristic of the large gap regions where the STM 
spectrum resembles that in the pseudogap phase d 12 ' 14 . 



Later, in the pseudogap phase, a similar checkerboard 
pattern was also observed 1 ^. Finally, the STM studies of 
Ca2-a;Na 2 Cu02Cl2 revealed the presence of a global long 
range commensurate checkerboard order independent of 
doping 1 ^. 

There have been various theoretical proposals to 
explain the non-dispersive checkerboard modulations. 
Most of these proposals are related to the competing or- 
ders. In these theories, the origin of the non-dispersive 
modulations are tied to the existence of particular order 
parameters. The theories including pair density mod- 
ulatior£ i 6 i 19 i 20 i 21 , current density modulatio n 22 ' 23 , spin 
modulation 2 , stripe charge modulation 7 -^, and impu- 
rity scattering 2 ^ and so on. 

Among the proposed mechanisms, the pair density 
wave (PDW) has been shown to capture important char- 
acteristics of the checkerboard density modulation. The 
mechanism of PDW derives from high pairing energy 
scale in cuprates. It suggests that, unlike the super- 
conductivity of normal BCS type superconductors that 
can be destroyed by breaking Cooper pairs, the super- 
conductivity in cuprates can be more easily weakened or 
destroyed by phase fluctuations than by pair breaking. 
Based on this argument, pair density localization* was 
first proposed to explain the local checkerboard modu- 
lation in the presence of impurity or vortex. Later, a 
global pair density wave(PDW) was proposed to explain 
the checkerboard physics in the pseudogap state 6 -^. It 
has also been shown that the symmetry of the tunnel- 
ing intensity can distinguish the pair density modulation 
from the conventional density modulation 6 . While the 
pair density modulation provides a good understanding 
of the experimental results, the theory does not cover 
two important characterizations of the tunneling experi- 
ments, namely, the dichotomy between nodal and antin- 
odal regions and the Fermic arc. 

Recently, we have proposed a d-wave checkerboard 
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density (DWCB) order in P-H channel 2 ^. The DWCB 
can be viewed as a natural extension of the d- 
density wave(DDW) order proposed to explain pseudo- 
gap physics^, and is only different from the latter in 
terms of order wavevectors. We have shown that the 
DWCB order must exist when the PDW order is present 
in the global d-wave superconducting state. Moreover, 
we have demonstrated that the DWCB captures many 
general features of the STM experimental results. It has 
been demonstrated that the DWCB order has little effect 
on the density of state at low energy in the superconduct- 
ing phase, but has a strong effect on the STM spectrum 
around the superconducting gap at high energy. This re- 
sult naturally explains the puzzling dichotomy between 
the nodal and antinodal regions observed in STM 14 and 
angle resolved photoemission spectroscopy (APRES) 27 . 
The DWCB order also preserves in FT-STS spectrum 
at the wavevectors, Q = {(7r/2a, 0), (0, 7r/2a)}, the same 
symmetry as that observed in the experiments. More- 
over, the DWCB preserves the nodes in the single particle 
spectrum, and generates a Fermi arc with little disper- 
sion around the nodal points at high temperature, which 
are consistent with the results from ARPES. The Fermi 
arc has been a signature of the pseudogap region, and 
has been proposed to explain the checkerboard pattern 
observed in the pseudogap stated. Thus, the DWCB 
provides a physical origin of the Fermi arc. 

In this paper, by solving the Bogoliubov-deGennes 
equation self-consistently, we further demonstrate the 
results drawn in ref.—. We design a microscopic 
model with the space modulated density-density inter- 
action. We show that from the self-consistent solutions 
of the Bogoliubov-deGennes equations, a weak spatially 
modulated density-density interactions can generate the 
DWCB and the PDW orders simultaneously in general. 
The combination of these two orders captures many im- 
portant features of the STM experimental results. We or- 
ganized the rest of this paper as follows: In Sec. II We in- 
troduce the d- wave checkerboard density (DWCB) order 
to study the electronic states in the disordered d-wave su- 
perconducting states, and the pseudogap phase at a tem- 
perature above T c . In Sec. Ill we show how the comple- 
mentary connection between the density orderings in the 
particle-particle and the particle-hole channels are closely 
related in the order wavevector and the symmetries. In 
Sec. IV we study the orders in cuprates by calculating a 
full self-consistent BdG equation including DWCB and 
PDW in the d-wave superconducting state. In Sec.V 
we calculate the local density of states and the spectral 
weights in the presence of both DWCB and PDW and 
compare the recent experimental data. 



II. d-WAVE CHECKERBOARD DENSITY 
ORDER IN CUPRATES 

In this secion, we shall introduce the DWCB order 
and discuss its roles in different regions of the phase di- 



agram in high temperature superconductors. In refi 2 ^, 
we have proposed a DWCB to explain the experimental 
results observed in the STM measurements in a global d- 
wave superconducting state. It has been shown that the 
presence of DWCB and competition with DSC can cap- 
ture the physics in a disordered superconducting state 
and in the pseudogap phase at a temperature above 
T c . Here we present a more detailed discussion and 
analysis. Firstly, we will show the average density of 
states(DOS) and the Fourier component at the order 
wavevector Q = {{ir/2, 0), (0, n/2)} calculated in a mixed 
state of DSC and DWCB orders show good agreements 
with experiments such as the spatially modulated LDOS 
in the disordered superconductors. Secondly, we will 
show the presence of the DWCB order with the d-wave 
superconducting order at a higher temperature above T c 
plays an important role in understanding the physics 
of the pseudogap phase: an emergence of the energy- 
independent Fermi arcs above T c and the linear depen- 
dence of its length on temperature. 

The mean-field Hamiltonian for the system where the 
DWCB coexists with d-wave superconducting order is 
given by 



Hmf — Ho+ Hdwcb j 



(1) 



where 



n a = &cLck CT + A k c k T c -ki + h - c - ( 2 ) 

k,cr k 

Hdwcb = Wk,q4 +qCT c k(T + h.c. (3) 

k,q,a" 

Tio is the mean-field Hamiltonian for the d-wave super- 
conducting state and thus Ak = Ao/2 (cosk y — cosk y ). 
The DWCB has particular order wavevectors, Q = 
{(7r/2,0), (0,7r/2)} and can be written as Wk, q = Wk<5 qi Q 
with Wk = Wo/ 2 (cosk x — cosk y ). 

To obtain a clear picture of the DWCB, we rewrite 
Hdwcb in the real space configuration: 



"Hdwcb = y^Im(Wp) 

r 

(cos — - s 
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where r = {x,y)a, Br iv) = J2A c l,a c r+ax( y ),a + h.c.) is 
the density operator defined in the nearest neighbor bond 
between r and r + ax(y) with x{y) unit vectors along 



x(y) directions, and J, 



x{y) 



r+ax(y) 



a — h.c.) is 



the current density operator defined in the same bond 
as Br^ y \ Fig. [1] shows a static pattern of the bond 
strength of the DWCB order, (B r ). It is clear that the 
DWCB order defined in Eq.Q has 4a x 4a periodicity 
and d x 2_ y 2 symmetry. Similar order parameters have 
been mentioned in Ref. 25 . 



FIG. 1: The configuration of the bond density of the DWCB 
order in the real space. It is manifestly shown that the pattern 
has 4a x 4a periodicity and d x 2_ y 2 symmetry. 

A. DWCB in the disordered Superconducting 
state 

Based on the mean- field Hamiltonian, Eq.JT]), we have 
calculated the averaged local density of states, p(ui), and 
the Fourier component of local density states(LDOS) at 
Q = {{tt/2, 0), (0, 7r/2)}, pq(uj), in two band dispersions. 
The calculation results turned out to be rather general 
and insensitive to the bare band structures. The effec- 
tive mean- field Hamiltonian, Eq.|T]), can be rewritten by 
using Nambu formalism: 

H MF =X>k#( k )V>k (5) 

k 

where ^ k = (c k t, Ck+QT, cL ki , c! k _ Qi ) , and 

/ 6c VF k A k \ 

U(V\ - W *- ^ k+Q Ak+Q (fi\ 
H W- A* -£_ k W£ +Q ( 6 ) 

V A* +Q Wic+Q -C-k+q/ 

Then the retarded Green function is given by 

G- 1 (k,iu) = (u J + tr 1 )I-H(k), (7) 

where I is the identity matrix with the same rank with 
H(k). The averaged density of states(DOS) and the 
Fourier component at Q can be calculated as the fol- 
lowing, respectively: 

p(w) = -iVlmG u (k, W ) (8) 

k 

Pq(u) = --VlmG 12 (k, W ) (9) 

k 

Firstly, we performed calculations in the particle-hole 
symmetric case, where the band dispersion is given by 

6c = — t/2 (cos k x + cos k y ) — p (10) 
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FIG. 2: (a) Averaged local density of states (LDOS) are 
plotted for various DWCB orders: Wo = 0, 4, and 8meV with 
t = — 125meV, t' = p = 0. (b) The Fourier components of 
LDOS at Q = {(tt/2, 0), (0, vr/2)}. 

We chose t — 125meV and p — 0. Ao = 40meV, which 
is relevant for underdoped BSCCO. The imaginary part 
of the self energy 77 = 5mcV was used for the numerical 
calculation. 

Fig. HJa) shows the averaged DOS normalized by one 
of the non-interacting Fermi liquids. In the absence of 
DWCB order, i.e. Wo = 0, there are sharp coherence 
peaks at the energy of superconducting gap, as expected. 
As DWCB order develops, the coherence peaks located at 
40mcV are suppressed, while the spectrum at low energy 
remains unchanged. From Wq = 4meV the prominent 
peak begins to appear within the superconducting gap. 
Note that even small DWCB made a strong effect on the 
spectrum at high energy as Wq increases. 

Fig. [U(b) shows the Fourier components of LDOS at 
the wavevectors, Q = {(tt/2, 0), (0, 7r/2)}. It is known 6 
that Pq(u>) = pq(-uj) for a bond-centered P-H pairing, 
while pq(u)) = — Pq(— uj) for a site-centered P-H pairing 
such as a conventional charge density wave(CDW). As 
expected from the fact that the DWCB order is bond- 
centered, Pq(u>) is even with respect to w, and it shows 
good agreement with experiment 9 . 

We repeat our calculations with the more realistic band 
dispersion provided by Norman et al^ and the result is 
displayed at Fig. [3] The band energy dispersion is now 
modified as such 

6c = ti/2(cos k x + cos k y ) + t<z cos k x cos k y 

+ f3/2(cos 2k x + cos 2k y ) 

+ i4/2(cos 2k x cos k y + cos k x cos 2k y ) 

+ t$ cos 2k x cos 2k y — p, (11) 
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FIG. 3: Averaged local density of states with the finite chem- 
ical potential included in the band dispersion provided by 
Norman et aS^. 



where t x = -0.5951eV, t 2 = 0.1636eV, t 3 = -0.0519eV, 
t 4 = -01117eV, and t 5 = 0.0510eV 29 . The chemical po- 
tential fi is now set to -0.1660cV. Compared with the 
P-H symmetric case, the effect of DWCB on the aver- 
aged DOS and the Fourier component of LDOS at Q is 
insensitive to the energy band structure. Qualitatively 
the numerical results are strikingly consistent with ex- 
perimental results^, and the large gap region can be in- 
terpreted in terms of the coexistence of weak(8-12meV) 
DWCB and DSC orders. 



B. DWCB in the Pseudogap state 

The results of LDOS have demonstrated the consis- 
tency between the coexistence of the DWCB and DSC 
orders and the STM experimental results in the disor- 
dered superconducting state. Now we will show that 
the presence of the DWCB order also captures important 
physics in the pseudogap phase at a temperature above 
T c . While there still have been hot debates over the inter- 
pretation of the origin of the the pseudogap phase, we will 
show how many experimental observations in this phase 
can be explained by interpreting the pseudogap phase as 
the mixed state of the DWCB and DSC orders. In this 
paper, we will focus on two features of the Fermi arcs 
developed from the nodal point along the Fermi surface: 
non-dispersive energy-independence and linear tempera- 
ture dependence. 

The autocorrelation of ARPES data from BSCCO 
show non-dispersive peaks in momentum space arising 
from the tips of the Fermi arcs in the pseudogap phase, 
while the superconducting state shows dispersion with 
binding energy-^. For the temperature dependence of the 
Fermi arcs in the pseudogap phase, it is known that its 
length increases linearly with temperature in the range 
between T c and T* , below which the material is believed 
to be in the pseudogap phase 30 . If the pseudogap phase is 
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FIG. 4: The spectral function given by Eq. (|12|l based on 
the energy dispersion, (|10|l . (a) The energy dependence 
of the Fermi arc along the Fermi surface. The white bars 
with the magnitude of |Q| = 7r/2 are displayed to show the 
Fermi arc is non-dispersive. The deviations for the energies 
{u> = 0, —10, —20, —30 meV) from |Q| are negligibly small by 
8%, 2%, 0%, and 13% of |Q|, respectively, (b) In the prsence 
of DSC(A = 40 meV) and DWCB(W = 8meV), the spectral 
weights, A{k.,uj = OmeV), are plotted with varying tempera- 
ture (T = 50K, WOK, 150K, and 20070 . 



strongly connected to phase fluctuations of ci-wave super- 
conductivity, the single particle spectrum should reflect 
the DWCB order. Therefore, a robust Fermi arc feature 
should exist in the mixed DWCB and DSC phases at high 
temperature. We will show that this is indeed the case. 

In order to illustrate the emergence of the non- 
dispersive Fermi arc in the pseudogap state with DWCB 
order coexistent with DSC, we have calculated the spec- 
tral function A(k, to) as given by the imaginary part of 
the retarded Green function: 



A(k,to) 



1 



ImGn(k, lo) 



(12) 



First we have studied the energy dependence of the Fermi 
arcs in the pseudogap state. We have plotted the spec- 
tral function A(k, lo) in the first Brillouin zone in Fig. [4] 
based on the energy dispersion, Eq. (fT0|) . In Fig. [4|a), 
the Fermi arcs are in red, and one of the d-wave checker- 
board wavevectors, Q = (7r/2,0), is shown as a white 
bar. The scattering wavevectors connecting the tips of 
each arc are nearly equal to the order wavevectors of 
DWCB, |Q| = 7r/2, which is consistent with the non- 
dispersive Fermi arc in experimental observations. The 
small dispersion is negligible when compared with the 
elongation of the gapless part along the Fermi surface 
from the nodal point in the DSC state without DWCB 
order. It may depend on the band width of the calcula- 
tion. When including more hopping terms in the energy 
band, the dispersion shown in the simple band calcula- 
tion will be reduced. 

Now let us consider the temperature dependence of 
A(k, lo) at the Fermi level(cj = OmeV). It can be calcu- 
lated in the pseudogap phase by taking the temperature 
dependence as an effect from the self energy, r\. We have 
plotted the spectral function A(k, lo) as a function of the 
temperature in the first Brillouin zone in Fig. H^b). At 
very low temperature the Fermi surface is gapped ex- 
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cept at the nodal point, (ir/2,ir/2). As the temperature 
rises, the nodal points grow significantly along the Fermi 
surface with slight broadening in the direction perpen- 
dicular to the Fermi surface. We will show the detailed 
calculations and discuss more in section V. 



tion, we will study how these density ordcrings in both 
channels are related in the modulation wavevector and 
the symmetry. 



The Connection between PDW and DDW 



III. COMPLEMENTARY CONNECTION 
BETWEEN THE ORDERS IN THE P-P AND P-H 
CHANNELS 

In this section, we will give a general argument regard- 
ing the complementary connection between the orders in 
the particle-particle and the particle-hole channels in the 
disordered (2-wave superconducting state at zero temper- 
ature as well as in the pseudogap phase at high temper- 
ature above T c . As shown above, a c?-wave checkerboard 
order (DWCB) in the P-H channel, (4 CT c k+Q(7 ) = $/(k) 
with Q = {(tt/2,0), (0,7r/2)} and /(k) = cos k x - cos k y , 
can explain the experiments on cuprates in the disordered 
DSC state and in the pseudogap phase. Preserving the 
same symmetry in FT-STS spectrum as that observed in 
experiments, the DWCB generates Fermi arcs with little 
dispersion with the binding energy around nodal points 
at high temperature above T c . Since the Fermi arc has 
been a signature of the pseudogap, the DWCB provides 
a direct link between the competing order and the pseu- 
dogap physics. 

The DWCB, however, is not a completely independent 
mechanism. In fact, there are some intimate connections 
with other orders in the P-P channel. Due to the exis- 
tence of a complementary connection between orders in 
both channels, the DWCB is directly connected to the 
PDW order. The angular symmetry of one orderings is 
determined by the combination of the angular symme- 
tries of DSC and the complementary order. In this see- 



To show the connection, we will use the example of the 
(i-density wave (DDW) order. Since the density order in 
the P-H channel has been well studied, examining the 
complementary connected ordering in the P-P channel 
will be a preliminary step to generalize the complemen- 
tary connection in the case of DWCB. 

The mean field Hamiltonian in the DSC state coexist- 
ing with DDW can be written as 

H = ^ £k4. CT Cka + iHkcJ^Ck+Qff + A k cj cT c^ ki + h.c. 



k<r 



= ^t (k)j4(k)V , (k)) 

where ^(k) = ( Ckt c k+QT c^ ki ct k _ Qi ) f , and 



(13) 



A(k) 



a iWk 

-iw k -a 

A* 
-AT 




(14) 



where c k(T and Ck<r are the creation operator and de- 
struction operator of the single particle with spin a, 
respectively. £k is the energy dispersion of the sin- 
gle particle. A k = A /2 (cos k x — cos k y ) and Wk = 
W /2 (cos k x cos ky) are the d-wave supeprconducting 
order (DSC) and the DDW, respectively. Then the 
Hamiltonian is diagonalized by the following transforma- 
tion: 



J 



/ CkT \ 

C k+QT 



u 



k-Ql 



-jWk 



£k+-Ek 



-iW k 



Jk-Sj, 



A /2B k (Bk-«k) ^2E k (E k +£ k ) ^/2E k (E k +^ k ) ^2E k (E k -(, k ) 

(, k -E k -iW k j k + E k -iW k 

y/2E k (E k -(, k ) y/2E k {E k +(, k ) y/2E k (E k +(, k ) y/2E k (E k -(, k ) 

^ Ak 

^/2E k (E k +^ k ) 

A k Q A k 

\ y/2E k (E k -(, k ) y/2E k (E k +e, k ) 



V / 2E k (E k -£ k ) 





/7i(k) \ 
72 (k) 
7a(k) 

V 74 (k) / 



r 



(15) 



where £ k = V Sk + l A k| 2 + W^. The eigenvalues E a (k) 
corresponding to the eigenvectors 7«(k) are given by 



tions: 



Si(k) =E 2 (k) = +£ k 
E 3 (k) = E 4 (k) = -E k . 



(16) 
(17) 



7 i(k)|0)= 72 (k)|n) = o, 

73 t (k)|^= 7 i(k)|O)=0. 



(18) 
(19) 



The ground state |fi) is defined by the following condi- 



Then non-zero expectation values of the bilinear opera- 
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tors of c's are 



( c kT C -kl) 



(4 T Ck+Qt) = 



|Akl 2 



\ c kT C -k-Qi/ 



where v 2 = \ [1 - jfc 
out DDW, or W k = 



Ak 
2E k 

2 iW k Ak 



(20) 
(21) 
(22) 
(23) 



k |A k p + iy k 2 
Note that in a DSC state with- 



= kT c k T) 



(24) 



Therefore, in the mixed state of the DSC, (c k T c — kl) oc 
Ak, and DDW, (c ko .c k+ Q CT ) °c iW k , a PDW order with 
the same wavevector Q is expected to exist naturally, 
whose symmetry is an extended s-wave manifestly given 
by 



( c kT c -k-Qi> K iM/ kA k 



iW A 



(cos k x — cos fe 9 



(25) 

This indicates that the mixed state of DDW and DSC 
leads to the existence of PDW. This argument can also be 
applied to the state of DSC coexisting with PDW leading 
to DDW. Therefore the mixed state of DSC and DDW is 
nothing but a complcmentarily connected description of 
the mixed state of DSC and PDW. 



B. 



The Connection between PDW and DWCB 



Given the close resemblance between DDW and 
DWCB, it is natural to expect that similar results also 
hold for the DWCB order, since the difference between 
DDW and DWCB is the order wavevector Q. Due to the 
high energy scale of superconducting gap in HTS, SC can 
be destroyed or suppressed by phase fluctuations. Car- 
rying out a similar analysis as Gossamer superconduc- 
tors^, we can show how the complementary connection 
between PDW and DWCB arises in the low energy ef- 
fective Hamiltonian that describes a system with pair 
modulation induced by phase fluctuations. 

We start from a BCS Hamiltonian on a 2D square lat- 
tice, 

Wbcs = ^2 £kC k!( r C k, CT + ^ Ak4,T c -k4 + hx -> ( 26 ) 

k,(r k 

The superconducting vacuum is constructed by the 
Cooper pairs with opposite momentum and spin: 

|*bcs) = «7bcs|0) oc e£* akC M o -*.i|0), (27) 



where |a k | 2 



and E k 



A local 



pair fluctuation induced by disorder such as a vortex, 
impurities, or other factors will lead to a new ground 
state with nonzero supercurrent due to the uncertainty 
principle. In other words, a new ground state \^f) can be 
obtained from the BCS ground state I^bcs) by applying 
a boost of the total momentum q of the Cooper pairs, 



I*) 



aV E k . q (*q4 T C -k + qi -*q C -k + qiC kT ) I 



BCS/ 



(28) 



where 77 is a small parameter and $ q is the structural fac- 
tor of pair fluctuations determined phenomenologically. 
For a commensurate checkerboard modulation, <i> q will 
be peaked at the related total momenta q of the pair, 
that is q = Q = {(tt/2, 0), (0, tt/2)}: $ q ~ <5 q , Q . 

The effective Hamiltonian associated with as a 
ground state can be obtained as, up to a first order 77, 



Hcfi = H.BCS + Hkii 



H 



hit 



(29) 



where Hkin and Hint are additional terms generated by 
boosted kinetic and pair interaction terms, respectively, 
in Eq. (|26p : for a commensurate checkerboard modula- 
tion, 



= Ek Ak, Q (c kiT c t _ k+qa + ct k+qjT 4 a ) + h.c, (30) 

where A k q = —i]J2 q £k&q ~ £k<S q ,Q and it represents 
an extended s-wave density order in the P-P channel, or 
PDW. And 

Hint = V X! *qAkC k ,o- c k+q,<x + h.c. 



E^4 



k, CT c k4 



h.c. 



(31) 



k,CT 



where Wk = r) <I> q A k ~ A k (5 qj Q and it represents a 
d-wave density order in the P-H channel, DWCB. The 
above derivation tells us that a e?-wave-like ordering in 
the P-H channel effectively leads to Cooper pairs with a 
finite center of mass momentum. 



IV. SELF-CONSISTENT CALCULATION FOR 
PDW AND DWCB 

Due to the complementary connection between PDW 
and DWCB, it is suggestive to calculate a fully self- 
consistent Bogoliubov-de Gennes (BdG) equation includ- 
ing both orders simultaneously. It has been argued that 
long-range interaction between charge carriers is very im- 
portant in cuprates. PDW can be induced by the long 
range interaction between hole-pairs. In this section, we 
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will show that a solution of a self-consistent calculation 
with modulating pair potential and nearest neighbor in- 
teraction can lead to orderings in both P-P and P-H chan- 
nels. 

We start from a full model Hamiltonian on a 2D lattice, 



(32) 



where m = + mi- 

In particular, we will consider nearest-neighbor hop- 
ping, tij = t and an attractive nearest-neighbor interac- 
tion, Vij. The chemical potential /i is chosen such that 
it is P-H symmetric case. The nearest-neighbor inter- 
actions between two opposite spins are considered to be 
attractive and modulating. They are given by a nega- 
tive constant pair potential and a long-range interaction 
modulated by Q = {(tt/2, 0), (0, tt/2)}, 



Vij = -Vo+SVy 

= -V + Y^ AV iJ ( COS Q • r i + C0S Q ' r J ) > ( 33 ) 
Q 



where 



+Vi, j = i±x 
-Vi, j = i±y 



(34) 



with Vq and V\ positive constants. While the local be- 
havior of SVij seems to be anisotropic at the center 
and (±2a, ±2a) as seen in Fig. El the global behavior 
is isotropic in a sense that there is no preferential direc- 
tion for attraction enhanced by negative modulation or 
weakened by positive modulation. The anisotropic mod- 
ulation is rotationally invariant with the combination of 
translation by 4a. 

Starting with this Hamiltonian, let us derive the BdG 
equations by introducing mean-field decoupling of the 
nearest neighbor interaction terms. The quartic term 
niaTija 1 has three different types of mean-field decou- 
plings 

C ia C i<j C ja' C ja' 
^ ( c icr C ia) C j<7' C jcr' ( C ja' C jcr') C i<7 C ia 

^~( C ia C jcr') C jcr' C ia C i<j C ia' ( C jcr' C ia) 



■0 



c ia c ja' i C ja' C ia 



(35) 



The first two terms are nothing but the chemical poten- 
tial terms. The decouplings in the second line are the 
usual BCS decouplings. The last two terms are decou- 
plings in the P-H channel, which are just the exchange 
terms. For interactions between two particles with same 
spins, these exchange terms effectively modified hopping 



+3V, 




FIG. 5: A modulating nearest-neighbor interaction, SVij, is 
plotted in a 2D square lattice. 



terms. For interactions between two particles with oppo- 
site spins, it corresponds to a bond-centered spin density 
wave. From our full Hamiltonian, we can argue that the 
spin rotational symmetry is not broken. Thus the SDW 
decouplings can be safely ignored. Also, since we are in- 
terested in d-wave pairing, the BCS decouplings of triplet 
pairing will be dropped, too. We then arrive at the fol- 
lowing MF Hamiltonian, 

HMF = [{t + Wij^Cjv + h.C. 



E [ A lf4 



(36) 



T c il + h - c - 



where 



Wi 
A 



(2) 



(i,3) 



-V ((CjiCit) + {CilCjl)) I 2 
-$Vij{c] a Cia) 

SVij ({cjiCit) + ( c il c jT>)/ 2 
H&ij - Vij{c\cj- a ) 



(37) 
(38) 
(39) 
(40) 



The pairing amplitude on a bond A ij\ stems from 

a constant nearest neighbor attraction, — Vq, and the den- 
sity wave pairings in the P-P and the P-H channels are 
caused by a long-range interaction, SVij. Since our major 
interest is to study the orderings in both channels at the 
same time, the chemical potential, Hij a , is taken to be a 
constant. 

The Bogoliubov-de Gcnnes equation is given by 



Ho A* 



-Ho 



u n (i) 
v n {i) 



= E n 



u n (i) 
v n (i) 



(41) 
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where Tto and A are transfer matrices such that 



n x(i) = -(E/* 

A -(*) = E, ( A S 



h A 



J„ T /XjiJ ( 42 ) 

?A <j) (43) 



where x(i) can be either u n (i) or v n (i). 

We numerically solve BdG eigenvalues E n and eigen- 
vectors (u n (i), v n (i)) on a lattice of N sites with a peri- 
odic boundary condition. We then calculate the pairing 
amplitude, and the density orderings in the P-P and the 
P-H channels which are, respectively give by, 

A ? = -yE^0X(0 + «»(iK(j)) ( 44 ) 

n 

A S } = yEWKW + ^M)) (45) 

n 

Wy = -*V«5>»0X(i). (46) 

Substituting these equations back into the BdG equa- 
tion, Ea.(|H]). and repeating the same process until self- 
consisitency is achieved for each of the local variables, 
we can obtain solutions of the BdG equations. In the 
calculation, we have used Broyden's method for efficient 
iteration. 

From the self-consistent solution, we can then com- 
pute the interesting order parameters. The cZ-wave pair- 
ing amplitude (DSC) at a lattice site r is determined by 
four nearest neighbors, Ay : 

A DSC (r) = ( A<J + . + A«_. - x <-> y)/A. (47) 

In the same way, an extended s-wave pairing density wave 
order in the P-P channels (PDW) can be obtained by 



A 



(2) 



Ap D w(r) = (a£> + A« + x «-> y)/4, (48) 



and cZ- wave density ordering in the P-H channel (DWCB) 
is 

WWcb(i-) = (lf r , r+i + W T ,r-s - x ^ y)/4. (49) 

We shall now present our results of full calculation. 
The uniform part of the interaction is set to Vb = — 2.5Z 
while the modulation part has amplitudes of V\ with the 
range of 0.1Z to l.Ot. The solutions we have obtained are 
independent of the initial guesses for the local variables. 
In Fig. [5] (a), the order parameters, AdscM, ApDw( r ), 
and Wdwcb(i"), in a 2D real space are plotted for Vq — 
2.5t and Vi = 0.5* with t = 125 meV. The P-H channel 
ordering, W T T > shows cZ-wave symmetry such that 



W T .r> = 



W (cos Q r + cos Q ■ r') , r' = r ± x 
-W (cosQ-r " ' 



where Wo = 0.06Z = 8.3 meV and the order wavevectors 
Q = {(±7r/2,0),(0,±7r/2)} as shown in Fig 1(b). 

The order in P-P channel has both d-wave, A™„ and 

(2) 

extended s-wave, A^,, part. The extended s-wave part 
is a pure modulation described by 

.(2) _ J Ai (cosQ ■ r + cosQ • r') , r'=r±x , > 
< Ai(cosQ-r + cosQ •!■'), r' =r±y [b ' 

where Ai = 0.1 3t = 16.3 meV and the same order 
wavevectors Q. We have also obtained the results of 



eZ-wave pairing amplitude, or, A 



i) 



Aq and A 



,1) 

r.r±y 



+ cos Q ■ r 7 ) , r 7 = r ± y ^ ' 



-Ao with Ao = 0.28t. In Fig. [7] the modulation ampli- 
tudes of A D sc(r), Apdw(i\), and Wdwcb(i") were plotted 
as a function of V% . 

The self-consistent BdG calculation shows that there is 
indeed a complementary connection between the density 
order in the P-P channel (PDW) and in the P-H channel 
(DWCB). Since PDW and DWCB are based on the same 
modulating interaction, SVij , not only PDW and DWCB 
should be present simultaneously with <i-wave supercon- 
ductivity but also the symmetries and order wavevector 
Q of both orders should be closely related to each other. 
Therefore, the mean-field solutions in the disordered d- 
wave superconductivity should include a cZ-wave density 
wave with Q in the P-H pairing and an extended s-wave 
density wave with the same Q in the P-H pairing. Fur- 
thermore it is important to take both density orders into 
account on the same footing when we try to understand 
the electronic states of the disordered high T c cuprates 
and the pseudogap state at high temperature. 



V. THE MIXED DWCB AND PDW STATE 

Analytically, the general features in STM measure- 
ments can be captured by the DWCB. Due to the 
anisotropy inherited from the (Z-wave factor of pairing, 
a weak DWCB order has a much stronger effect on the 
antinodal region than on the nodal region. Thus it natu- 
rally explains the puzzling dichotomy between the nodal 
and antinodal excitations in high temperature supercon- 
ductors. The local phase fluctuations of Cooper pairs 
lead to a local modulation of cZ-wave ordering in the P-H 
channel (DWCB), which strongly affects the antinodal 
single particle excitations, as well as an extended s-wave 
order in the P-P channel (PDW). Since both PDW and 
DWCB are bond-centered, Pq(u>) will be an even func- 
tion of uj 6 . This symmetry distinguishes the PDW and 
DWCB orders from the typical CDW in the P-H channel. 
The existing experimental results are consistent with the 
even case. 

The purpose of this section is to determine the symme- 
tries of density orders in P-P and P-H channels and then 
study the pseudogap physics under the presence of the 
orderings in P-P and P-H channels. Due to the comple- 
mentary connection, if the PDW is an extended s-wave, 
DWCB is a <i-wave P-H pairing. If the PDW is a d-wave, 
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FIG. 6: (a) The amplitude of the order parameters, Adsc(V), Apdw(t")i an d Wdwcb(t") defined by Eq. (|47l) . (|48|l . and ((49}, 
respectively, plotted in 16 x 16 lattice sites. Adsc(?") ~ Ao, Apdw = Ai cos Q r, and W / DwcB(r) = Wo cos Q -r with Ao = 34.2 
meV, Ai = 16.3 meV, and Wo = 8.3 meV, respectively, t = 125 meV is chosen to have the amplitudes in units of meV. (b) 
The Fourier transforms of the order parameters are displayed in the first Brillouin zone. The prominent peaks are located at 
Q = {(±7r/2, 0), (0, ±tt/2)} in both PDW and DWCB. 




8.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 7: The amplitudes, Ao, Ai, and Wo of A D sc(r), 
ApDw(r), and WdwcbM, respectively, as a function of the 
modulation long range interaction Vi. t = 125 meV is chosen 
to have the amplitudes in units of meV. For a d-wave super- 
conducting order (DSC), Ao, denoted by x along the curve 
for DSC, at each Vi represents the mean value of AdscM 
which is slightly modulated. 



the P-H pairing density order must be extended s-wave. 
It is worth noting that the symmetries of the orders are 
very important to explain the behavior of the conduc- 
tance spectra at low energy. We will show that in DSC 
state the d-wave P-H pairing(DWCB) and the extended 
s-wave P-P pairing(PDW) can explain the low energy 
spectra while the solution with other symmetries cannot 
sufficiently match the experimental results. 



In section II, we have shown that, even without PDW, 
the presence of DWCB and DSC orders in a pseudogap 
state can capture the features of the Fermi arc. To un- 
derstand the temperature dependence of the length of 
the Fermi arc above T c we assumed the temperature de- 
pendence of the imaginary part of the energy 7/ in the 
Green function implicitly. Using the results of the self- 
consistent calculations, this assumption of the tempera- 
ture dependence of 77 is not necessary. The presence of 
the complementary connected c£-wave orders in P-P and 
the extended s-wave order in P-H channels with DSC or- 
der naturally lead us to the temperature dependence of 
the length of the Fermi arc within the Franz and Millis 
models. 
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FIG. 8: (a) The averaged DOS are plotted for A D sc = 33 
meV in the case of d-wave P-H and extended s-wave P-P 
channel orders: A PDW = 12.7 meV, W DWC b = 3.3 meV. The 
arrow indicates Sh, that is, DOS at u — off of one for pure d- 
wave superconducting state. Sh = 0.0599. (b) The averaged 
DOS in the mixed state of DSC (Adsc = 33 meV) with d- 
wave P-P and extended s-wave P-H channel orders. Here we 
plotted one with A PDW = 17.3 meV, and W SWCB = 10 meV. 
Sh = 0.2933. 



A. The symmetries in P-P and P-H channels 

In this subsection we study how the local density of 
states is associated with the symmetries of the comple- 
mentary connected orders in P-P and P-H channels in a 
disordered d-wave superconducting state. 

The local density of states(LDOS) in the DSC state in 
the presence of P-P and P-H channel orders is calculated 
by 

p(r,w) = J2[\ u n(r)\ 2 S(u-E n ) + \vn(r)\ 2 S(u+E n )} (52) 

n 

The averaged DOS is given by 

p(u) = ^2p{T,u) (53) 

r 

and the Fourier components at the wavevectors Q = 
{(7r/2,0),(Q,7r/2)}by 



PQ 



M=]T e ^Xr,c) 



(54) 



We have calculated them in a simple band structure, t = 
— 125meV, and p = 0. 

For the d-wave P-H and extended s-wave P-P pairing, 
we have used Vq = 2.5i and various V\ as in the ta- 
ble ([55]) . where order parameters as solutions of the BdG 



calculation are given in meV, and Sh is the height of p(u) 
at u> = off of the superconducting DOS displayed in a 
dashed line in Fig. [5J 



Vl 


Adsc 


Apow 


Wdwcb 


Sh 


0.4* 


33.0 


12.7 


3.3 


0.0599 


0.6* 


35.7 


19.9 


9.9 


0.1125 


0.8* 


39.6 


27.5 


13.0 


0.1445 


1.0* 


44.2 


35.6 


16.1 


0.1546 



(55) 



For the extended s-wave P-H and d-wave P-P channel 
orders, we have used Vq = 2.5* and V\ as in (|56p. where 
W^SWCB denotes the extended s-wave P-H pairing for con- 
venience. 





Adsc 


ApDW 


WSWCB 


Sh 


0.4* 


31.5 


11.5 


6.6 


0.1144 


0.6* 


32.7 


17.3 


10.0 


0.2933 


0.8* 


34.2 


22.9 


13.3 


0.4200 


1.0* 


36.2 


28.4 


16.6 


0.3589 



(56) 



The BdG calculation shows that DWCB, which is d- 
wave P-H pairing order, is crucial to the physics in the 
disordered superconductor. It captures the experimental 
observation of conductance spectrum at both low and 
high energy. As expected, the d-wave symmetry of the 
P-H pairing density order has little effect at low energy 
while it strongly affects the superconducting coherence 
peaks. Therefore we can conclude the competing order 
in the P-H channel must have d-wave symmetry. Due to 
the complementary connection in the orders in P-H and 
P-P channels, it naturally lead to the existence of the 
pair density order that has extended s-wave symmetry. 

In Fig.[8j to illustrate the above argument, we plotted 
the averaged DOS in two cases: (a) d-wave P-H and ex- 
tended s-wave P-P channel orders, and (b) s-wave P-H 
and extended d-wave P-P channel orders in the d-wave 
superconducting state (Adsc = 33meV). 



B. The DWCB and PDW in the pseudogap state 

Now let us consider the spectral function A(k, w) in the 
pseudogap state to capture the features in the Fermi arc. 
First we rewrite the mean-field Hamiltonian for PDW 
and DWCB in momentum space. Taking the Fourier 
transform of Eq. (|36[) , we obtain 



HpdW = X] A ri+s4,*4+6,-<T + h " c - 



(57) 



E 



(2) t 



-i(p+q)-r e -iq-<5 



t t 
p,o- u q,-CTi 



where S denotes the unit vectors for the nearest-neighbors 
of r. Since A^ +(5 = Ai/2 £q[ cos Q • r + cos Q ■ (r + $)], 
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the term in the square bracket in Eq. (I57|) is rewritten as 
Ai t cos Q • r + cos Q ■ (r + S^e-^+^e^ 5 

r,Q,<5 
Q,<5 

= Ai y^[cosp a + cos(p x + Q x ) + x <-> y]S p+cl .Q 
Q 

Thus Hpdw m momentum space can be given by 

Wpdw = £ Ai 2) 4 )<r ct k _ Q) _ ff + h.c. (58) 

k,Q,er 

where A J, = Ai cos /c^ -|- cos (k x + Q x ) +x<-yy\. 

The momentum space expressions for DWCB in 7Ymf, 
Eq. (|36j) . can be obtained in the same way: 



k,Q,er 

where Wk = Wo [cos fec+cos^a; +Q X )- 
the mean-field Hamiltonian, Eq. ([36 
momentum space as 



h.c. 



(59) 



x <-» y] . Therefore 
, can be recast in 



Hmf = 



E 

k,cr 



A (1) c f c 



EG 

k,Q,<r 



(2). 
'k 



,t 

k.rT 1 - 



-k-Q,- CT + Wwcl +Qa Ck t , 



and A 



(i) 



where £k = — t/2(cosk x + cosk y ) — fx, auu ^ k 
Ao/2(cosfc a; — cosky). In the same way as in section II, 
it can be re-expressed in terms of the Nambu formalism 
as such 



where ipu 



k 

(ckT,c k +QT,cL ki ,cL k _ Qi ) , and 



(61) 



/ a W k A k 1} A k 

W k * Ck +Q A k 2 ] Q 

.(!)* A (2), f 

^k ^k+Q S-k 

\A k 2) * A k % W k+Q C-k+q/ 



,(2) 



JJ(k) = 



\ 



^k* + Q 



(62) 



The spectral function is given by the imaginary part 
of the Green function, 



Gii(k,w) = (fil^ + ^-^k)]^ 1 ^) 



(63) 



The ground state |f2) is defined by the state without 
quasiparticles. In the pseudogap phase, it can be cal- 
culated within the model of Franz and Millis 32 : 



A(k,w) = -i J du'P(uj')ImG n {k,u}-u'), 



(64) 




(60) 
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FIG. 9: The energy dependence of the the Femi arc. (a) 
shows the development of the gapless region as a function of 
energy u> in the superconducting state without DWCB. As 
u increases (from top to bottom), the scattering wavevectors 
connecting one tip of the green curve to another tip decrease. 
However, in (b) the gapless region remains unchanged as oj 
increases, that is, the Fermi arc is non-dispersive. 



where oj' is the energy shift due to the presence of a 
uniform supercurrent, and P(u>') is the probability dis- 
tribution of oj' calculated within a 2D XY model as such 



P(cj') = V2nW(T)e 



-u' 2 /2W 2 (T) 



(65) 



where W(T) is an order of pseudogap and an increasing 
function of temperature above T c . 

Now we will study the energy dependence of the Fermi 
arc in the pseudogap state. One of the salient features 
of the pseudogap phase is that the Fermi arc is non- 
dispersive, whereas it is dispersive in the superconducting 
phase. The Fermi arc along the Fermi surface is identi- 
fied by the peak of EDC at oj = 0. Fig. [^b) demon- 
strates that the Fermi arc in green remains unchanged 
as oj is increased. In the superconducting state, however, 
the gapless region develops as the energy decreases below 
the Fermi level as seen in Fig. [^a). 

While the Fermi arc is non-dispersive, it is found that 
the length of it is increasing linearly as a function of tem- 
perature 3 ^. Within the model of Franz and Millis, when 
T > T c , W(T) is an approximately linearly increasing 
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FIG. 10: The temperature dependence of the Femi arc above T c . (a) EDC's along the Fermi surface from the nodal point (top) 
to anti-nodal point (bottom) are plotted. The arrows indicate the gapless part in the Fermi surface and show the length of the 
Fermi arc linearly depends on the temperature in the pseudogap state, (b) The Fermi arc as a function of the temperature T 
at above T c . 



function when the order parameters remain fixed. Adopt- 
ing the values of W(T) above T c , we have calculated 
A(k, u>) as a function of temperature. 



In Fig. fTUf a). EDC are plotted from the nodal to antin- 
odal point along the Fermi surface. Clearly it shows a 
gapless nodal point and gapped anti-nodal point as in- 
dicated by darker curves. But, contrary to a supercon- 
ducting state where the Fermi surface is gapped except 
at the nodal point, there is a region, the so called Fermi 
arc, where peaks still survive at u> = and thus the 
superconducting gap is closed. The curve on top is the 
EDC at the nodal point, (7r/2, 7r/2), perpendicular to the 
Fermi surface. It is manifest that the size of the arrows 
is increasing as the temperature rises. Fig. [TOlb) shows 
the Fermi arc in the first Brillouin zone. As in the Fig. [3] 
the green color represents the gapless region of the Fermi 
surface, that is, the Fermi arc. 



VI. CONCLUSION AND DISCUSSION 



We have presented a detailed analysis of competing or- 
ders in cuprates which explains the checkerboard pattern 
observed in STM spectra. We have shown that in general 
there is a complementary connection between the orders 
in the P-P and the P-H channels in both the supercon- 
ducting state and the pseudogap state. The symmetries 
of the orders in both channels are closely related to each 
other. In the d-wave superconducting state, the pres- 
ence of the d-wave(extended s-wave) P-H density order 
implies the coexistence of the extended s-wave(ci-wave) 
P-P pairing order. The self-consistent calculations in the 
disordered superconducting state result in solutions with 
uniform d-wave superconducting order, a P-H channel 
checkerboard density order with the order wavevector Q, 
and a P-P channel pairing density wave with the same 
Q. The symmetries of the channel orders is crucial to 
the effects on the density of states in STM experiments. 
We have found that d-wave checkerboard density order 
in the P-P and the extended s-wave pair density wave 
are in good agreement with the experimental data in the 
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superconducting phase. 

The coexistence of orders in both P-P and P-H chan- 
nels are important to understand the pseudogap phase 
which can be considered as a state which maintains pair- 
ing amplitude without phase coherence. The fluctuation 
of phase would naturally lead to the DWCB and PDW 
orders. We have shown that the presence of DWCB and 
PDW results in a non-dispersive Fermi arc and our cal- 
culations show the linear dependence of the length of the 
Fermi arc at temperatures above T c . Moreover, the ef- 
fect of the presence of PDW and DWCB on single particle 
spectra is much larger in the antinodal direction than in 
the nodal direction. These results also explain the di- 
chotomy between the nodal and antinodal single particle 
spectra in the cuprates. 



The connections of the P-P and P-H channel orders are 
important in formulating the effective low energy theory 
in cuprates. From our results in this paper, the low en- 
ergy effective theory should include both orders. So far, 
most theories have treated them independently. It is also 
clear that the presence of both orders can result in new 
physics in transport and thermal properties. A detailed 
study of these effects will be reported elsewhere. 
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